INHOMOGENEOUS AND COUPLED CONTINUOUS TIME 
RANDOM WALK LIMITS 



PETER STRAKA 



Abstract. We study stochastic process limits of Continuous Time Random Walks 
(CTRWs) whose jumps and waiting times vary in both time and space and are possibly 
coupled. We prove limit theorems for the weak convergence in Skorokhod space of 
cadlag paths. In the case of coupled waiting times and jumps, there exist two different 
limit processes. These limit processes are then characterized via a tuple of coefficients, 
which is derived from a jump-diffusion processes in space-time. We then calculate the 
one-dimensional laws of these two limit processes, and show that these laws solve 
forward and backward Kolmogorov equations with non-local time operators. Finally, 
we discuss an application of our theory to the case of subdiffusion in a time-dependent 
force-field, anomalous diffusion and functional limit theorem and fractional derivative 
and subordination and coupled random walks and fractional kinetics 



I. Introduction 

Continuous Time Random Walks (CTRWs) are generalizations of random walks 
which allow for the modelling of long resting periods, infinite variance jumps, and the 
coupling thereof. CTRWs are hence versatile models for anomalous diffusion processes 
and have found applications in a variety of fields such as physics [1] , biology [2] , finance 
[3] and hydrology [I]. There is a rich interplay between scaling limits of CTRWs and 
fractional differential equations: Waiting times whose tail probability has parameter 
a G (0, 1) usually lead to diffusion equations in which the time derivative has order 
a [S], and to the time-change by the inverse of a stable subordinator. For the theory 
on fractional differential equations of this kind see e.g. [HI El IH] - If the jumps have a 
regularly varying probability distribution with parameter /3 e (0,2), this yields frac- 
tional derivatives in space [9] . The functional limit theory for stochastic process limits 
of CTRWs has been developed in the series of papers [TQl [HI EE21 [13] for the cases of 
uncoupled and coupled waiting times, and governing differential equations were derived. 

This paper aims at further generalizing the limit theory to coupled and uncoupled 
CTRWs for which the distribution of jumps and waiting times depends on space and 
time. Such processes are of interest in the physics literature [HI [H)] , and arise naturally 
e.g. in fractional chemotaxis problems [IB]: Say particles of two species A and B are 
moving randomly according to subdiffusion, and A is attracted by B. This is modelled 
via a drift of A along the concentration gradient of species B, which changes in time 
and space. The "standard approach" of changing the clock process of A to an inverse 
stable subordinator then fails, even in the uncoupled case; this is particularly evident 
if the external time-dependence is non-random [17] and has in the past led to some 
confusion in the physics literature j!8j . 
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A theory for inhomogeneous jumps and waiting times has been developed in [19], with 
a focus on the position-dependence of jumps and waiting times. We extend the theory of 
functional limit theorems, limit laws and governing differential equations further, with a 
focus on time- dependence and a coupling between jumps and waiting times. In section 
[21 we prove a functional limit theorem for scaling limits of CTRWs. As discovered 
independently in [20] and [15] . there are two different limit processes (CTRW limit 
and "overshooting" CTRW limit) if long waiting times occur simultaneously with long 
jumps. We also show how the limit processes can be represented as jump diffusions in 
space-time, yielding a convenient representation via a tuple of coefficients. In section [3] 
we calculate the laws of the two limit processes using the theory of entrance- and 
exit-distributions of Markov processes. In section H] we derive backward differential 
equations for the two limit processes, and achieve in some sense generalized versions of 
Kolmogorov's backward equations. In section 0, we derive forward differential equations 
for the case where the coefficients of the jumps depend on both time and space. We 
close with an application of the developed theory to a subdiffusion process in a space- 
and time- varying potential. 

2. Stochastic Process Limits of CTRWs 

Continuous Time Random Walks (CTRWs) in state space M d are given by 

i) a sequence J\, J%, . . . of M. d - valued random variables called the jumps and 

ii) a sequence W\, W2, ■ ■ ■ of (0, oo)-valued random variables called the waiting 
times. 

We consider scaling limits of CTRWs, and hence we introduce a scaling parameter 
t > 0, which we will attach to every random quantity. We then define the processes 

[t/rj [tM 

(2.1) A t = H J ii D i = E w *> 

k=l k=l 

which have the interpretation that A T {nr) is the position of the walker after n steps, 
and D T {nr) is the sum of n waiting times. Now consider the process inverse to D\, 
which is defined as 

(2.2) El = inf{s > : D T (s) > t} 

and is also called the hitting time process. It takes values in rN, and if El = rk then 
the number of waiting times which have elapsed at or before time t is k — 1. If one hence 
assumes that the k-th jump JJ, of the walker occurs at the start of the k-th waiting 
time Wjl, then this yields the process 

(2.3) Y; = A T (El), 

and similarly if JJ, occurs at the end of the k-th waiting time then resulting process is 

(2.4) Xl = A\El-r). 

In the literature, XI is commonly called a Continuous Time Random Walk (CTRW) 
and dates back to Montroll and Weiss [21]. The process Y t T has been of interest in the 
recent literature [221 [13], an d has been termed overshooting CTRW (OCTRW). 
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The scaling limit theory for XI and Y t T has been studied for i.i.d. jumps and waiting 
times p3] and in the general case in [20]. We consider weak limits of probability 
measures on Skorohod path space D(R d ), that is the set of all cadlag paths in W 1 , 
endowed with Skorohod's J-topology [23] . 

In the setup as above, the limit theorem reads as follows: 

Theorem 2.1 ([20]). Let DJ be non-decreasing and unbounded for all r > 0, and 
assume that as r j. the K d+1 -valued process (Al,D^) converges weakly on (B(IR d+1 ), J) 
to a limiting process (A t , D t ), where D t is unbounded and has strictly increasing sample 
paths a.s. Then the process E t inverse to D t has continuous sample paths a.s., and 
the processes X[ and Y t T converge weakly on (B(IR d ), J) to limiting processes X t and Y t 
given by 

(2.5) X(t-) = A(E t —) Y t = A(E t ). 

(Note: Both X t and Y t are cadlag, and in general A(E t —) ^ A(E(t—)) = Y(t—).) 

We will refer to X t and Y t by CTRWL and OCTRWL. The two limit processes X t 
and Y t are equal e.g. in the uncoupled case, where each waiting time is independent of 
each jump. In the coupled case, however, X t and Y t are different processes in general. 

In this paper, we study CTRWs and OCTRWs which are renewed after every pair 
(J£, W£) of jumps and waiting times. To be precise, we assume that the discrete process 

(2-6) k^(Al Tl Dl T ) 

is a Markov chain in space-time M. d+1 whose law is defined by a an initial distribution 
and a jump kernel K T (x, z; dy, dw) such that 

(2.7) P[j; +1 G dy, W T n+l G dw\A T nT = x, D T nT = z] = K T (x, z; dy, dw). 

The paths of CTRWs and OCTRWs are uniquely determined by each realization of 
the above Markov chain, and we have the following interpretation: The CTRW X[ is 
renewed after every jump, and the OCTRW Y t T is renewed after every waiting time. 
The space-time process (A t ,D t ). As r decreases to 0, we assume that the jumps and 
waiting times become smaller, resulting in smaller and more frequent jumps. We in- 
troduce a class of potential limit processes (A t , D t ) as follows: Let b : M. d+1 — > M. d , 
7 : M. d+1 — > M + and a : M, d+l — > M. dxd be bounded continuous functions, where a(z,x) 
takes values in the non-negative definite matrices. We interpret these as coefficients 
named "spatial drift" , "temporal drift" and "diffusivity." We also define a "space-time 
jump kernel" K(x,z;dy,dw) as a kernel from M. d+1 to {[0,oo) x M. d } \ {(0,0)} which 
satisfies 

(2.8) J K{x,t;dy,dw) ((||y|| 2 + \w\) A l) < oo, (x,t)eR d+1 , 

where A denotes the minimum of two entities. We assume that (x, t) h->- J g(y, w)K(x, t; dy, dw) 
is continuous whenever g(y, w) is a bounded continuous function which equals in a 
neighborhood of the origin. Now let {T s } s >o be a Feller semigroup, acting on the Ba- 
nach space Co(M. d+1 ) of real- valued functions defined on which are vanishing at 
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oo, and assume that T s admits an infinitesimal generator of the form 



Cf(x, t) = bi(x, t)dif(x, t) + j(x, t)d t f(x, t) + - ^2 ayj(x t t)d iS f{x, t) 

(2.9) 



+ 



f(x + y,t + w)- f{x, t)-^2 hi(y)dif(x, t) 



i=l 



K(x, t; dy, dw) 



where the vector (x, t) G R d+1 is indexed by 1, . . . , d, d+1 and where hi(y) = yil{\\y\\ < 
1} is a truncation function. We then assume a stochastic basis (Q, J 7 , jF s , F x,t ) which sat- 
isfies the usual conditions (i.e. right-continuity of J 7 , and augmentation by P x, *-null sets) 
and an adapted process (A,, D s ) on this basis whose finite-dimensional distributions are 
governed by T s and the Chapman-Kolmogorov equations. In particular, the measure 
¥ x,t only charges sample paths starting from (x,t) (that is, F x,t ((A , D ) = (x,t)) = 1), 
and we denote expectation with respect to F x,t by the operator ~K x,t . Processes satisfying 
such conditions exist, see e.g. [2U Ch.6]. 

The following proposition gives sufficient conditions for the convergence of {A T S , D T S ) 
to (A S ,D S ): 

Proposition 2.2. Suppose the following convergence statements hold locally uniformly 
in (x, t): 

(2.10) limr" 1 / h i {y)K T {x,t;dy,dw) = b i {x,t), l<i<d 

TiO J 

(2.11) limlimr -1 / wK T {x,t;dy,dw) = ~j{x,t) 

eiO r^O J 

{(y,w):0<w<e} 

(2.12) limliniT -1 / yiyjK T (x,t;dy,dw) = aij(x,t), 1 < i,j < d 

eiO tIO J 

{{y,w):\\y\\<e} 

(2.13) limr" 1 / g(y,w)K T (x,t;dy,dw)= / g(y,w)K(x,t; dy, dw), 

J J 



for each real-valued bounded continuous function g(y, w) defined on M. d+1 which equals 
in a neighborhood of the origin (0,0). Assume moreover a family {?7 t }t>o of probability 
measures on M. d+l such that lim r ^ ? 7 T = V holds in the weak sense. Then (Al,D"[) 
with initial distribution rf converges weakly on (D(M d+1 ), J) to the process (A t , D t ) with 
initial distribution rj, as r \. 0. 

Proof. We apply theorem IX. 4. 8]. First, we replace time t in {A\,Dl) with the 
process rN(t/r), where N t is an independent standard Poisson process of intensity 
1. Observe that this has no effect on X t T and Y t T : the trajectories of the CTRW and 
OCTRW only depend on the set of points traversed by {A\, DJ) and not on the time at 
which these points are traversed. By slight abuse of notation, we again write (Aj, DJ) 
for the process resulting from this time-change. We have thus achieved that (AJ, DJ) 
is a pure jump Markov process in continuous time, with state space Its jumps 
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arrive at the rate r 1 . Relative to the truncation function 

fx, w) if II fx, w)\\ < 1 



/i(x, w) 



else 



the process (Al,DJ) is hence a semimartingale with characteristics {{B T , T T ), C T , v 1 
relative to h, where 



t 

-l 



Bl it) = J bl(A T s , D T s )ds, 6[ (x, «;) = r" 1 y fcfo, v) K T (x, w; dy, dv) 

T T (t) = J* 7 r (^I, D T a )ds, Y(x, w) = r- 1 J h d+1 (y, v) K T {x, w; dy, dv) 

c ij(t) = J c ij( A l, DT s )ds, 4,(x, w) = r' 1 J hihj{y, v) K T (x, w; dy, dv) 

u T (ds; dy, dv) = K T (A T S , D T S ; dy, dv)ds. 
Observe that (12.111) and (I2.12p are equivalent to 

(2.14) limr -1 / wK T (x, z; dy, dw) = j(x, z) + / h d+1 (y,v)K(x,w; dy, dv) 
T ±° J J 

{(y,w):0<w<l} 

(2.15) 

limr -1 / yiyjK T (x, z;dy,dw) = a,ij(x, z) + / hJij(y,v)K(x,w;dy,dv), 
J J 
{{y,w)-\\v\\< 1 } 

1 < i,j < d, an d that relative to h, the d + 1-component of the first characteristic of 
(A t , D t ) is 

(2.16) T t = J -/{A, D s )ds + J h d+1 (y, v)K(A„ D s ; dy, dv), 

and hence our requirements match with the those of the cited result. □ 



Together with theorem 12 .1\ the above proposition shows that a large class of stochastic 
process limits of CTRWs and OCTRWs of renewal type can be effectively expressed via a 
tuple of coefficient^ (7, a, b, K). Moreover, proposition ^. 2l provides sufficient conditions 
on a sequence of (O) CTRWs to be an approximating sequence of an (O)CTRW limit 
process with given coefficients (7, a, b, K). 

3. Laws of CTRW Limits 

In this section we calculate the laws of the CTRW limit X t and the OCTRW limit 
Y t for fixed t. The main idea is to formulate our setting as an entrance-exit problem 
for a Markov process: X t is the last position of A before (A s , D s ) exits from the set 
M. d x (—00, t], and Y t is first position of A when (A, D) enters the set R d x [t, 00). 



^We note however that this representation is subject to one degree of freedom: for c > 0, the tuple 
(7, a, b, K) will yield the same CTRW and OCTRW limit as (07, ca, cb, cK). 
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We introduce the O-potential kernel U of (A,D), which arises from the semigroup 
(T s ) and acts on the non-negative measurable functions / defined on ~R d+1 as follows: 

(3.1) 

/roo |~ poo 

f(y,w)U(x,t;dy,dw) = / T s f(x, t)ds — E x '* / f(A(s),D(s))ds 
Jo Uo 

Note that U f can be infinite, and hence we make the assumption 
(3.2) 

(x,t) i — y Uf(x,t) is bounded if / is supported on R d x [a, b], — oo < a < b < oo. 

This assumption reflects the "transience" of D s , and is satisfied e.g. if D s is a subordi- 
nator [26]. Next, we introduce the random set C of points t G K which are traversed 
continuously by D s : 

C := {(t, u) G K x : t = D r (w) = A--M for some r > 0} 

We write {£ G C} or C t for the event {u G f2 : (£, w) G C}, and we calculate the law of 
X t and Y t separately on the two sets {t ^ C} and {t G C}. 

Proposition 3.1. Fix (x, s) G IR d+1 and let t > s. Then 



E*'-[/(X t _,y t );{^C}]= y U(x,s-dy,dv) J f(y,y + z)K(y,v;dz,dw) 

R d x[s,t] R d x[t-v,oo) 

for all bounded measurable f defined on M. d xM. d . 

Proof. The complement of C(u) in R is a.s. a countable union of closed intervals 
[D U -,D U ], where u is a jump instant of the process D. We have Ae 4 - < t < -Ds t 
for all t > s. One also checks that AD^ t > on the set {t ^ C}. On the set {t G C}, 
one has ADs t = 0: To see this, first note that {E t < u} = {D u > t} a.s. and hence E t 
is an ^-stopping time. On the set {t G C}, the sample path of E t is (left-) increasing 
at t, and hence the J-"-optional time E t is announced by the optional times E t _i/ n . E t 
is hence J-"-predictable [27J p.410]. Without loss of generality we choose the canonical 
setting Q = D(lR d+1 ) with natural filtration T t . Then {A t ,D t ) is quasi left-continuous 
[271 prop.22.20], and A(A,D) Et = a.s. 

Thus writing J = {(r, u) G R + x Q : A(A, £>) r 7^ (0, 0)} for the optional set of jump 
times, we find 

f(X t .,Y t )l{t i C} = £)/(4-, A r )l{D r _ <t<D r } 

(3.3) reJ 

= f( A r-, A- + AA r )l{D r _ < t < A_ + AD r }. 

The last expression equals J W(cu, r; y, w)/i(uj, dr; dy, dw) for the optional random mea- 
sure \x on R + x R d+1 associated with the jumps of the process (A, D), 

(3.4) fi(u,dr;dy,dw) = ^ 1 j( v ^) S (v,A(A,D) v (u 1 ))(dr;dy,dw), 

v>0 



and the predictable integrand 

W(u,r;y,w) 

:= f(A r _{u), A r _(u) + y)l{D r _(uj) < t < D r _{u) + w}. 
The compensator ji p of fi equals 

(3.5) /U p (w; dr, dy, dw) = K(A r _{oj), D r _{uj); dy, dw)dr. 

According to the compensation formula [251 H.1.8], taking expectations on both sides 
of ( 13. 3p then yields the proposition. □ 

Proposition 3.2. Fix (x,s) G Suppose that the measure U(x,s;-) is Lebesgue- 

absolutely continuous with density (y,w) H- u(x,s;y,w). Then for Lebesgue- almost 
every t > s and every bounded continuous f , 



(3.6) E x >'[f(Y t )l{teC}]= / f(yHy,t)u(s,x;y,t)dx. 

Jm. d 

Moreover, the above holds true ifY t is replaced by X t -, X t or Y t -. 

Proof. Similarly as in [28], D r admits a decomposition into a continuous and a discon- 
tinuous part via 

D c r = ( j(A v ,D v )dv, D d r = £ AD V , r>0. 

To see this, we first note that (A r , D r ) is a semimartingale, and hence D r allows the 
decomposition 

(3.7) D r = ^AD V 1{AD V > 1} + B r + M r 

v<r 

where B r is a predictable process of finite variation (the first characteristic of D r ) and 
M r is a local martingale. Due to [251 IX §4a] and ( 12.91) . B r = J Q r j(A v , D v )dv where 
j(y,w) = 7(y,w) + J vl{\\(z,v)\\ < l}K(y,w;dz,dv). Since D r has no diffusive part, 
M r is purely discontinuous and equal to 

M r = J2^D v l{AD v < 1}- II iwl{||(2/, <l}K{A v ,D v -dy,dw)ds. 

v<r 

But then ([331) reads D r = D d r + D c r . 

For fixed u, the paths of D., D c and D d are non- decreasing and hence define Lebesgue- 
Stieltjes measures dD, dD c and dD d on [0,oo). Then for any bounded continuous / 
and g we have 

poo /*oo 

(3.8) / f(A v )g(DMA v ,D v )dv= / f(A v )g(D v )dD c v . 
Jo Jo 

The continuous measure dD c does not charge the countable set {v : AD V ^ 0} of 
discontinuities of D. and coincides with dD on the complement {v : AD V = 0}. Hence 

7 



the last integral can be written as 

POO 

/ f{A v )g{D v )l{v: AD v = 0}dD v . 
Jo 

Using the substitution formula 

poo poo 

/ f(r)dD r = / f(E r )dr, 
Jo Jo 

valid for all strictly increasing (not necessarily continuous) D and continuous /, we get 

POO 

/ f(Y t )g(D Et )l{t : AD Et = 0}dt. 
Jo 

Now note that AD Et = is equivalent to t £ C and implies D Et —t. Hence the above 
lines show that the left hand side of (13.81) equals 

POO POO 

/ f(A v )g(D v )j{A v , D v )dv = / f(Y t )g(t)l{t £ C}dt. 
Jo Jo 

Applying expectations and Tonelli's theorem yields 



f(y)g(t)j(y,t)u(x,s;y,t)dydt = / W' a [f{Y t )l{t £ C}]g(t)dt. 

R d +! Jo 

Since g is an arbitrary bounded continuous function, this yields ( 13. 61) for almost every 
t. For any fixed t, we have already noted in the proof of the preceding proposition that 
A(A,D) Et = (0,0), P x ' ,s -almost surely on {t £ C}. This and the continuity of sample 
paths of E t shows X t -X t _=0 = Y t - Y t _ on {t £ C}. □ 

We can now characterize the joint law of X t and Y t at a fixed time t: 

Theorem 3.3. Fix a starting point x £ M. d and a starting time s £ R. Suppose the 
^-potential U(x,s;-) has a Lebesgue density (y,v) \-> u(x, s;y,v) which is continuous 
on the set M d x (s, oo). Then for every t > s and every bounded continuous f defined 
on R d x R d , 



(3.9) 

E*> s [f(X t , Y t )] = [ f{y, y)j{y, t)u(x, s; y, t)dy 



+ / / u(x,s;y,v) I I f(y,y + z)K(y,v;dz,dw)dydv. 



' w>t—v 

Moreover, Y t has no fixed discontinuities, and if K(y,v;M. d x {w}) = for all y,v,w, 
then Xt has no fixed discontinuities. 

Proof. Equation (13.91) holds for almost all t > s, which follows from the addition of 
Propositions 13.11 and I3.2| note that the domain w >t — v may be replaced by w > t — v 
since K(y, v; M. d x {w}) > for only countably many w > 0. Consider now a sequence 
t n — > t such that t n > t and (13.91) holds for t n . Then pointwise in Lebesgue-almost every 
y (resp. (y,v)), the integrands converge. Equation (13. 9p hence follows for all t > s, by 
right-continuity and the dominated convergence theorem. 



For the remaining statements, let 



G := {(t,u) 6lxH: D~(u) = t < D r (u) for some r > 0} 

denote the random set of left endpoints of the intervals contiguous to the range {D r : 
r > 0} of D, and let 

D:={(t,w)elxll: D~(u) < t = D r (u) for some r > 0} 

denote the random set of right endpoints. Recall that the process E has continuous 
non-decreasing sample paths a.s., since D is strictly increasing a.s. These sample paths 
are right- increasing and left-constant at t whenever t G D, and they are left-increasing 
and right-constant at t whenever t G G. If t > s is a discontiuity of the process X, i.e. 
{t G disc(X)}, we find A~(E(t)) ^ A~(E(t+)), which is equivalent to E(t) e disc(A) 
and E(t) being right- increasing at t. We have seen in the proof of Proposition 13. II that 
E(t) ^ disc(A) a.s. on the set {t G C}. Hence disc(X) C D a.s. Similarly, one sees 
that disc(Y) C G a.s. Now an application of the compensation formula yields 



F x ' s (t e D) = E 1 



H D r < t}i{AD r = t - d;} 



r>0 



E 2 



' poo 

/ 1{D; < t}K (A;, D~; R d x{t- D~}) dr 
Jo 



which vanishes under the conditions of Statement 2, and similarly we find 



F x ' s (3r >0:D: 



t<D r -e) 



E 1 



H D 7 = t}l{AD r > e} 



r>0 



1{D~ = t}K (A;, D~; M. d x (e, oo)) dr 







since D is a.s. strictly increasing. But since e > is arbitrary, this means F XyS (t G G) = 
0. □ 

We now define the transition kernels corresponding to the processes X and Y as 

(3.10) P s J(x) := E*''\f(X t )], Q s , t f(x) := E*>'[f{Y t )]. 

We can give simplified formulas for probability densities of X t and Y t : 

Corollary 3.4. Assume the conditions of Theorem \3.3[ Then the kernel P s ,t(x, dy) has 
a density P s ,t(x, y) whenever t > s, given by 



(3.11) P Sjt (x,y) = j(y,t)u(x,s;y,t) + / u(x, s; y, v)K(y, v; R d x (t - v, oo)) dv. 

For the kernel Q s ,t{x, dy), we find 
(3.12) 

Q s ,t(x, dy) = j(y, t)u(x, s; y, t)dy + 



u(x, s; z, v)K(z, v; (dy — z) x (t — v, oo)) dz dv. 



v=s 
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The above formulas extend results which were given in (T2J [13] for the case of constant 
coefficients and 7 = 0. 

4. Backward differential equations 

In this section we derive backward differential equations for the transition proba- 
bilities P s>t (x,y) and Q s ,t{%,dy) of CTRW and OCTRW limit processes. We interpret 
(x, s) as the backward and (y, t) as the forward variables. 

Theorem 4.1. Suppose the ^-potential kernel U(x,s;-,-) of the space-time process 
(A r ,D r ) admits a Lebesgue density (y,t) i-> u(x,s;y,t) for every (x,s) G M d+1 , which 
is such that the map (x,s) H- u(x,s;y,t) has bounded partial derivatives up to second 
order outside every neighborhood of (y,t). Moreover, assume that 

(4.1) C X:S J J u(x,s;y,t)F 1 (y,t)dydt = J J £ x>s u(x,s;y,t)F 1 (y,t)dydt, 

(4.2) £ XyS I u(x,s;y,t)F 2 (y)dy= / £ x , s u(x, s;y,t)F 2 (y) dy 



for all locally integrable functions F\ and F 2 for which the integrals are defined. Then 
for fixed t G R, 

(4.3) -£ x ,,Ps, t f{x) = f(x)K(x,s;R d x (t-s,oo)), s < t, P t , t f{x) = f(x), 



(4.4) -C x , s Q s ,tf(x) = J f(x + z)K(x, s; dz x (t - s, 00)), s<t, Q t , t f(x) = f(x), 
for every x E M. d and f G C (lR d ). 

Proof. We will apply the operator H x s to P Sit f(x) = J f(y)P S)t (x,y)dy. Our first step 
is to note that 

(4.5) £x, s u(x, s; y,t) = 

for Lebesgue almost all (y, t) G R d+1 and all (x, s) ^ (y, t). To see this, let G C7 (M d+1 ) 
have compact support. Then Ucp defines a bounded function due to assumption (13.21) . 
For A > 0, the A-potential of the semigroup T r is defined as the kernel / 1— > U x f = 
joo e -\r Tr f dr It jg known j29] that U x maps ^(M^" 1 " 1 ) into the domain of C, and that 

(A - £)U x <p = <p. But then 

CU<b = £ lim U x 6 = lim CU X 6 = lim XU X 6 - 6 = -6, 

XiO A4.0 A|0 

using the fact that £ is a closed operator. If the support of <f> does not contain (x, s), 

then 

(4.6) 



J J £ x>s u(x,s;y,t)(f)(y,t)dydt = C XjS J J u(x,s;y,t)<f)(y,t) dydt = £U(f)(x,s) 
Our next step is to write 



[x,s 



P s ,tf(x) = h(x, s) + J 2 (x, s) + I 3 (x, s), 
Qs,tf(x) = h{x, s) + J 4 (x, s) + I 5 (x, s), 
10 



where 

h(x,s)= I f(y)^(y,t)u(x,s;y,t)dy 



I 2 (x,s) = / f(y) / u(x, s;y,v)K(y,v;R x (t - v, oo))g(v) dvdy 

J y J v 

I 3 (x,s)= / f(y) / u(x, s;y,v)K(y,v;M d x (t-v, oo))(l - g(v)) dvdy 



h(x,s)= g(v)u(x,s;y,v) / f(y + z)K(y,v;dzxdw)dydv 

J y J v J z J w>t—v 



h(x,s)= / (1 - g(v))u(x,s;y,v) / / f(y + z)K(y 1 v;dzxdw)dydv 

J y J v<t J z J w>t—v 

and where g is a non-increasing continuous function defined on IR such that = 
(—00, s + e] and 5 ,_1 ({0}) = [i — e, 00), and where s<s + £:<t — e < t. Note that 
u(x, s;y,t) is not supported on {s > t}, hence we can omit the integration boundary 
v > s. By our assumption on u(x, s; y, t), we find Ch(x, s) = CI 3 (x, s) = CI^x, s) = 0. 
We can express 

(4.7) I 2 (x, s) = UF P (x, s), h(x, s) = UF Q (x, s) 

where F P (y, v) = f(y)g(v)K(y, v; R d x (t - v, 00)) and F Q (y, v) = g(v) J z J w>t _ v f(y + 
z)K(y,v; dz x dw) define functions in Co(M ). Hence we find 

CI 2 (x, s) = CUF P {x, s) = -F P {x, s) = -f(x)K(x, s; R d x (t - s, 00)), 
£/ 4 (x, s) = CUFq(x, s) = —Fq(x, s) = — J f(x + z)K(x, s; dz x (t — s, 00)). 

This shows the first statements in (j4.3p & f!4.4p . and P t t f(x) = f(x) and Qt,tf(%) = f(%) 
follow from the right-continuity of the sample paths of (A r ,D r ) and from D r having 
strictly increasing sample paths. 

□ 

Remark 4.2. The assumption that the potential measure of the space-time process 
(A r , D r ) has a twice continuously differentiable density with bounded derivatives seems 
to be a bearable restriction, but it is surprisingly hard to check. In [30J, condi- 
tions are given under which the transition probability of (A r , D r ) has a joint density 
p r (x, s; y, t) of class C 2 ; however, in order to transfer this property to the potential den- 
sity u(x, s;y,t) via a dominated convergence argument, one needs the uniform bound- 
edness of these derivatives, and boundedness is not shown in this book. 



5. Forward differential equations 

In this section we derive governing differential equations in the forward variables (y, t) 
for the (O)CTRW limit densities P s ,t(x, y) and Q Syt {x,dy). We rely on two simplifying 
assumptions: 
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(1) The law of the waiting times only depends on the current position in space x, and 
not the current time t. That is, we can write K(y, v; dz x dw) = K(y; dz x dw), 
and j(y, v) = j(y). 

(2) Jumps and waiting times are uncoupled. That is, the measures dz x dw t— > 
K(x; dz x dw) are supported only on the "coordinate axes" {R d x {0}) U ({0} x 
(0, oo)). This means that CTRW and OCTRW are indistinguishable, and 

P Sjt (x,y)dy = Q Stt (x,dy). 

The generator C in (I2.9p then decomposes as C = A + V. The operator A acts only 
on the space variable, but its coefficients may depend on time t, which we highlight 
by the notation A(t). Similarly, T> acts only on the time variable, but its coefficients 
may depend on space, and we write T>(x). If it is unclear which variable an operator is 
acting on (i.e. y oy z), we will put this variable in lower script, e.g. A y {t). The operators 
A and V, explicitly, are 



^ d 1 ^ 0^ 

A(t)f(x) = J2k(x,t)—f(x) + a^x,t)-^-f(x) 

(5.1) 4=1 1 *' i=1 1 ' 

r d 



+ 



f(x + y)- f(x) - hi(y)dif(x 

i=l 



K{x,t;dyx{0}), 



(5.2) V(x)g(t) = j(x)^-g(t) 4 [g(t + w) - g(t)]K(x; {0} x dw). 



dr 



w>0 



Hence if / G dom(A(t)) for all t and g G dom(T>(x)) for all x, then the function 
/ ® g(x,t) := f(x)g(t) lies in dom(C) and 

(5.3) jr,[f®g](x,t) = (A(t)f(x))g(t) + (V(x)g(t))f(x). 

Finally, we introduce a memory kernel Ai as follows: For every y G 7BL d , write dt i— > 
I(y, dt) for the measure on (0, oo) whose Laplace transform is 

poo poo 

(5.4) / e- xt I(y,dt) = \ 1 (y) + \h(y,\), h(y,X) = e~ xt K (y; R d x (t, oo)) dt. 
Jo Jo 

Such a measure exists; it is the 0-potential of the subordinator with drift ^(y) and Levy 
measure with tail function t y K(y; M. d x (t, oo)), see e.g. [26]. The operator A4. then 
acts on the transition probability P s .t{x, y) via 

(5.5) [MP] s , t (x,y) = P s , t - r (x,y)I(y,dr). 

Note that we extend P s> t{x,y) by the value whenever t < s, and hence the integral 
may be taken over the entire domain r G R. 

Theorem 5.1. Assume the conditions from Theorem \3.3\ and the beginning of this 
section. Then the transition densities P s j{x,y) satisfy the forward equation 

(5.6) P s J(x) = f(x), §- t P ^f^) = [MP] s ,tA(t)f(x), t > s, 
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whenever f £ f] t€R dom(A(t)) C Co(K d ) and x £ R d are snc/i £/iat t i— >■ P st f(x) has a 
continuous derivative for t > s and (x,t) > A(t)f(x) is uniformly continuous on ~R d+1 . 

Proof. Let g £ dom(V(x)) for all x, and let g have compact support contained in (s, oo). 
Then f ® g E dom(£), and 

f®g = U x {\-£){f®g) = \U\f ®g)- U x £(f ® o), 

where A > and t/ A denotes the resolvent operator associated with £ [3TJ. Letting 
A i together with (13.21) then shows that —U£(f ® g) exists and equals / <E> g. We 
then find 

-f®g = U£(f ®g) = U[(Af)g + (Vg)f], 

(5.7) -f(x)g(s) = (u(x, s; y, T),A v (T)f(y)g(r))y lT + (u(x, s; y, T),V T (y)g(r)f(y)) y>T , 

where we have used the notation (f,g)y, T = ff f(y,r)g(y,r)dydT. For the first term 
on the right-hand side of (15.71) . we use the formula 

(5.8) [MP] s , T {x,y) =u(x,s;y,r). 

To see that it holds, we first find the Laplace transform of Pq, t (x, y), 

POO 

P 0>T (x, y)e~ Xr dr = [ 7 (y) + h(y, A)] / u(x, 0; y, r)e" Ar dr, 
o Jo 

which shows 

POO PT 

Po,T-r(x,y)I(y,dr)= / u(x,0;y,r)dr, 



after an application of the convolution formula for Laplace transforms. A similar cal- 
culation with replaced by s then shows (15. 8p . For the second term on the right-hand 
side of (15. 7p . we consider a test function <p £ C£°((s, oo)), and calculate 

(5.9) 

(u(x, s; y, r),V(f>(r)) T = (u{x, s; y, r), 7 ( 2 /)0'(r) + jf [0(r + w) - (f>{r)]K{y; R d x dw)^ 



u(x, s; y, r), 7 (y)0'(r) + / <j>\r + w)K(y\ R d x (w, oo)) 

Jo 

u(x, s; y, r)y(y), <f>'(r)) T + / J u(x, s;y,r- w)K(y; R d x (w, oo)), 0'(r) 
P s>T {x, y), <p'(r)) T = - (JI~P s>T (x, y), 0(r)^> . 



Now we replace g in (15. 7p by a sequence of functions which approximate the Dirac 
distribution concentrated at t > s. Using (I5.8P and ( 15.90 . we arrive at 

(5.10) = ([MP] s , t (x, y),A y (t)f(y)) y - U^P.,t(x, y), f(y) 
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in the limit. This shows the sought differential equation, and P S}S f(x) = f(x) follows 
from the right-continuity of sample paths of (A r ,D r ) and D r being strictly increasing 
a.s. □ 



The above theorem extends the main result in [T3], which was the derivation of a 
Fokker-Planck equation in the case where D r is a stable subordinator and A r a process 
with continuous sample paths. It applies to space-dependent waiting times and a limit 
process A r which may have discontinuous sample paths. 

In the case where D r has i.i.d. increments (i.e. a Levy process started at s) and the 
jumps are homogeneous in time (i.e. A(t) does not depend on t), we have a simplified 
result (compare [191 Eq.(35)]). Define the formal adjoint of T>(x) as the operator 

(5.11) V\x)g(t) = --y( x )g{t) + / [g(t - w) - g(t)}K(x; {0} x dw). 

at Jw>0 

It satisfies 

(5.12) (V t (x)h(x, t), </>{x, t)) Xtt = {h(x, t),V*(x)(j>{x, t)) Xyt 
for every test function G C£°. 

Corollary 5.2. Assume the same conditions as in the above theorem, and additionally 
suppose that A r is homogeneous in time and D r is homogeneous in both time and space; 
that is, j(x,i) = 7 = const, K(x, t; {0} x (it;, 00)) = h{w) and A{t) = A for some 
operator A and all t. Then the transition densities P s ,t{x, y) of the CTRW limit satisfy 

(5.13) P s J(x) = /(x), -V* t P s J(x) = P s ,tAf(x) + f(x)h(t -s), t>8 
for each f G dom(A) and x EW 1 such that 1 1— > P s>t f(x) is differentiable on (s, 00). 
Proof. Without loss of generality let s = 0. Integrating (I5.6P with respect to t yields 

(5.14) P 0tt f{x) - f{x) = j P 0tt _ r Af(x)I(dr). 

Now let g(t) be an integrable function which is continuous on [0, 00) and which vanishes 
at 0. If its Laplace transform equals A 1— > g(X), then the Laplace transform of —T>*g(t) 
equals A (7 + h(X))g(\). Applying this, the definition (15.41) of I(dr) and uniqueness of 
Laplace transforms to the above equation yields the corollary. □ 

A similar result has been given for fractional dynamics, corresponding to a stable 
subordinator D r , see [Sj. 



6. Examples 

6.1. Subdiffusion in a space- and time-dependent force field. We consider a 
CTRW limit process that results from i.i.d. waiting times with Pareto distribution 
of tail parameter a G (0,1). Let ip(w) = a(l + u;)" 1- ", and introduce the scaling 
parameter A > 0. The rescaled density ip\{w) := tp{w/X)/X satisfies 

(6.1) A-"Va(*) ^ Mt) -= a*" 1 "" 
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As in [14], we consider nearest-neighbor jumps on a lattice with spacing Ax > 0. With 
probabilities r(t,x) and £(t,x), the particle at x jumps right or left at time t. These 
probabilities are given by a space- and time dependent potential V(x,t) via 

r(x, t) = C' 1 exp{-/3[V(x + Ax, t) - V(x 7 1)]}, 
£(x, t) = C~ l exp{-/3[V(x - Ax, t) - V(x, t)]} 

where C is defined by £(x,t) + r(x,t) = 1. If after k waiting times and jumps 
(Ji, Wi), . . . , (Jfe, Wk) the particle is at (x,t) in space-time, then the distribution of 
the k + 1-st waiting time and jump pair (Jk+i, Wk+i) is given by 

F((W k+1 , J fe+ i) edw x dy) = K T (x,t;y,w) dydw, 

K T (x, t; y, w) = if)x(w)[£(x, t + w)S(y + Ax) + r(x, t + w)S(y — Ax)] 

since the k + 1st jump occurs at time t + W^+i- The above can be simplified by noting 
that a process with the same distribution results if the k + 1-st jump is assumed to 
occur at time t and the starting time s is replaced by s + W , where W has the same 
distribution as all Wi. Now the kernel K T reads 

K T (x,t;y,w) = ip x (w)[£(x } t)8(y + Ax) + r(x, t)5(y - Ax)]. 

We choose the following dependence of the temporal scale A and the spatial scale Ax 



on r: 



(6.2) 



A' 



T 



r i 



Ax 



2 2 

~ a t 



a 



where T(-) denotes the Gamma function and a 2 is the diffusion constant and where 
/ ~ g means that the quotient of / and g tends to 1. 

We now check equations ( 12.101) to ( 12.131) to find that proposition 12.21 is fulfilled, with 

0, b(x, t) = -f3a 2 d x V(x, t), a(x, t) = a 2 , K(x, t; y, w) = 5(y)w^ a /T(-a). 
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Since now K(x,t;-,-) is supported by the coordinate axes, the corresponding CTRW 
limit and OCTRW limit are the same. The backward equation for the density (x, s) H- 
Ps,tf(x) then reads 

(6.3) 



d a 



2 dx 2 



Pa 2 



d_ 

dx 



P 8 , t f{x) = f{x 



(t 



r i 



a 



P t>t f(x) = f(x) 



where d a /d(—s) denotes the Riemann-Liouville negative fractional derivative of order 
a [19]. For the forward equation, we first note that D is a subordinator with symbol 
^(A) = X a (26] and the Laplace transform of I(y,dt) = I(dt) is \~ a . Inversion yields 
I(dt) = t Q_1 c?t/r(l — a). Hence convolution with the measure I is equivalent to a 
fractional integral of order a, and the memory kernel acts as a Riemann-Liouville 
fractional derivative of order 1 — a. The forward equation for P s>t (x,y) hence reads 



(6.4) ^Ps,tf{x) 







l-Q 



Of 



1-a 



P. 



a 2 d 2 | 2 d d_ 
2 dx 2 dx ' dx 



/Or), P 8t ,f(x) = f(x). 
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6.2. Space-dependent Levy walks. Levy walks were introduced in [32], and have 
only recently been studied on the stochastic process level [33] . Their main feature is that 
via a strong coupling between waiting time and jump, a finite travel "velocity" of the 
random walker is achieved at each step. Note however that each jump is instantaneous, 
and the term velocity is only to be interpreted in an averaged sense. The theory 
developed in this paper allows for a stochastic framework in which the velocity of the 
walker may vary in space. 

We introduce the space-dependent travel speed c(x) > 0, which we assume smooth 
and bounded away from both and oo. We specify the CTRW limit process via its 
parameters 

(6.5) a = 0, 6 = 0, 7 = 0, K(x; dz x dw) = ~r^N ( x dw 

c(x) \c(x) 

where N(dz x dw) is the Levy measure 

(6.6) N(dz xdw)= / 8 w6 (dz)l{w > 0}—^-\{d6) 

J S d-l W L + a 

on M d+1 . Here 5 y denotes the Dirac measure concentrated at y £ M d , and A is a spectral 
measure on the unit sphere S^ 1 in M. d . A Levy process in M. d+1 with characteristics 
(0, 0, N) is tightly coupled, in the sense that each jump in the M d -component is ac- 
companied by a jump in the M-component which has the same Euclidean length. In 
transitioning from N(dz x dw) to K(x; dz x dw) the waiting times are left unchanged, 
while the length of the spatial jumps is adjusted by the factor c(x). 

A sequence of CTRWs which converges to the limit process as defined above is 
achieved through letting r j 0, where for each r we consider a CTRW correspond- 
ing (for instance) to the space-time kernel 

f (x dw 

(6.7) K T {x; dz x dw) = 5 wc(x)e (dz)l{w > Oj -^jj^ - wr - 1/a y +Q A W 

6.3. Waiting times with space-dependent heavy tail parameter. In [31], CTRWs 
with varying "anomalous exponent" were studied. We assume that the heavy tail pa- 
rameter a(x) is a smooth function of space, which takes values in the interval (0, 1) and 
is bounded away from both and 1. Jumps are of nearest neighbor type and symmetric, 
on a lattice with spacing r 2 . Both waiting times and jumps are homogeneous in space, 
i.e. parameters do not depend on t. The space-time random walk then has transition 
kernel 

(6.8) K T (x- dz x dw) = -{5 T {dz) + 5^ T {dz))-_ 



2 V rv ' TK "Y{\- a(x))T 1 /«(*)(l + 1WT - 1 /«(*))i+«(*)' 
and in the limit r — > this yields the coefficients 

oc ( x i dw 

(6.9) a = l, 6 = 0, 7 = 0, K(x; dz x dw) = — 5 (dz). 

1 (1 — a{x))w l+a[ - x > 

In [31], a discrete lattice was considered on the bounded interval (0,1), with reflect- 
ing boundary conditions. It has been found that in the long-time limit the CTRW 
limit process is localized at the lattice point % at which a(x) attains its minimum, a 
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phenomenon called "anomalous aggregation" . The above provides a continuum setting 
for the work in [31], and it is an interesting open problem to confirm that localiza- 
tion persists in continuous space. The backward and forward governing equations for 
the CTRW limit are then the fractional differential equations of variable order. The 
backward equation is written as 

(6.10) — — -^P s t f(x) = —P s t f(x) + f(x)±- '—— -, P t t f(x) = fix). 

The memory kernel Ai in forward equation is symbolically written as the fractional 
derivative of variable order c 



Q t l-a(x) ■ 
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